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Figure  1.   Station  plot  of  MIZPAC  77  cruise.   The  position         3 
of  the  ice  margin  at  the  time  of  observation  is 
also  shown. 

Figure  2.   Distribution  and  intensity  of  f i nestructure  during       5 
MIZPAC  77.   Symbols  are  described  in  the  text. 
Isotherms  (°C)  are   the  maximum  temperature  in  the 
water  column. 

Figure  3.   Station  80  salinity  and  temperature.   The  solid        10 
line  was  obtained  while  lowering,  the  dashed  line 
while  hoisting.   The  dotted  line  is  the  result  of 
despiking  the  down  trace  but  not  smoothing.   See 
Figure  k   for  the  smoothed  curve. 

Figure  h.      Station  80.   The  down-going  salinity  profile  is        14 
on  the  right  as  a  solid  curve.   The  dashed  curve 
is  the  result  of  despiking  and  smoothing.   The  two 
curves  to  the  left  give  the  contributions  of  the 
short  time  constant  and  long  time  constant  correc- 
tions to  the  despiking.   See  Figure  3  for  the  un- 
smoothed  corrected  curve  and  to  see  the  temperature 
transients  which  caused  the  salinity  spikes. 


THE  OCEANOGRAPHIC  CRUISE  OF  USCGC  BURTON  ISLAND 
TO  THE  MARGINAL  SEA- ICE  ZONE  OF  THE  CHUKCHI  SEA- 
MI  ZPAC  77 

by 

Robert  G.  Paquette  and  Robert  H.  Bourke 
Naval  Postgraduate  School,  Monterey,  CA 


I.   INTRODUCTION 

This  report  presents  the  data  and  briefly  describes  the  oceanographic 
results  of  the  cruise  of  USCGC  BURTON  ISLAND  into  the  region  of  the  sea- 
ice  margin  of  the  Chukchi  Sea  during  the  period  2k   July  to  6  August  1977 
as  part  of  the  program  designated  MIZPAC  77.   The  primary  objective  of  the 
cruise  was  to  find  and  characterize  f i nestructure  in  the  vertical  tempera- 
ture profiles  and  to  discover  its  horizontal  distribution  and  causes.   This 
is  the  fifth  cruise  devoted  to  this  general  problem.   Other  cruises  in  1971, 
1972,  and  1 97^  were  reported  by  Paquette  and  Bourke  (1973,  1976)  and  MIZPAC 
75  by  Zuberbuhler  and  Roeder  (1976). 

I  I .   GENERAL  DESCRIPTION 

The  scientific  group  boarded  BURTON  ISLAND  at  Nome,  Alaska  on  the 
afternoon  of  2k   July,  one  day  later  than  in  initial  planning.   The  scientists 
and  their  affiliations  were: 

Dr.  Robert  G.  Paquette,  Naval  Postgraduate  School,  Chief  Scientist 

Dr.  Robert  H.  Bourke,  Naval  Postgraduate  School 

LCDR  Gordon  P.  Graham,  Canadian  Forces,  Student  at  Naval  Postgraduate 

School 
Mr.  Jonathan  D.  Trent,  Naval  Postgraduate  School 
Mr.  Peter  Becker,  Applied  Physics  Laboratory,  University  of  Washington, 

Seattle. 

The  scientific  group  disembarked  at  NARL,  Barrow,  Alaska  on  the  afternoon 
of  6  August. 

The  measurements  made  were  salinity  and  temperature  profiles  through- 
out the  entire  water  column  at  157  stations,  using  the  Applied  Physics  Labora- 
tory portable,  hand  lowered  CTD.   The  lowering  rate  of  the  CTD  was  about  1  m 
sec  "^  resulting  in  a  data  rate  of  approximately  three  points  per  meter.   The 
latter  was  checked  systematically  with  Nansen  bottles  lowered  on  a  second 
wire.   Prior  to  leaving  each  station,  the  temperature  and  salinity  were 
plotted  utilizing  a  Hewlett-Packard  9100  series  computer/plotter  system. 
These  rough  plots  were  used  to  make  immediate  assessments  of  the  presence 
of  fi nestructure  and  to  aid  in  the  decision  of  where  to  make  the  next  few 
stations. 


A  current  meter,  which  was  Intended  to  orient  lines  of  closely  spaced 
stations  along  the  flow  direction,  was  found  to  be  of  little  utility  be- 
cause of  the  effects  of  the  ship's  iron  on  the  magnetic  compass.   However, 
the  failure  to  orient  sampling  lines  along  a  presumed  direction  of  propa- 
gation was  at  least  partially  overcome  by  running  east-west  lines  as  well 
as  north-south  lines  in  areas  containing  f inestructure. 

Because  of  the  loss  of  one  day  of  cruise  time,  no  measurements  were 
made  in  Bering  Strait  and  the  southern  Chukchi  Sea  until  the  latitude  of 
Pt.  Hope  was  reached.   Otherwise,  we  proceeded  as  in  the  general  plan,  ex- 
ploring for  f inestructure  along  the  ice  margin  from  70°N  toward  Barrow  and 
studying  it  both  in  longitudinal  and  lateral  sections  when  found.   This  re- 
sulted in  much  intensive  study  near  Barrow  followed,  near  the  end  of  the 
cruise,  by  exploration  near  71°N  and,  finally,  a  few  more  measurements  near 
Barrow. 

III.  DATA 

The  CTD  was  standardized  by  means  of  a  Nansen  bottle  lowered  on  a 
second  wire  to  a  depth  just  above  the  sea  floor.   Forty-nine  such  compara- 
sons  were  in  sufficiently  unchanging  water  for  temperature  standardization 
and  kO   for  salinity.   The  mean  CTD  temperature  and  salinity  errors  were 
0.08°C  and  -0.021 °/00,    respectively.   The  standard  errors  of  these  means 
were  0.0051°C  and  0.005°/oo  and  the  standard  deviations  were  0.036°C  and 
0.025°/oo.   A  correction  was  applied  to  the  depth  data  to  account  for  the 
difference  in  density  between  fresh  and  salt  water,  as  the  CTD  pressure 
sensor  was  calibrated  in  fresh  water. 

The  CTD  records  its  data  on  a  cassette  which  is  eventually  transferred 
to  a  seven-track  tape  by  APL-UW  for  data  editing  and  analysis  at  NPS.   A 
computerized  editing  routine  was  written  to  remove  erroneous  data,  inter- 
polate data  where  necessary,  make  temperature,  salinity  and  depth  corrections, 
and  remove  spurious  salinity  spikes.   The  despiking  and  editing  procedures 
are  described  in  some  detail  in  Appendix  A. 

Heading  data  for  each  station  are  listed  in  Appendix  C.   These  contain 
station  position  and  number,  date/time  of  CTD  lowering,  water  depth,  type 
of  navigation,  wind,  wave,  and  air  temperature  data,  etc.   Appendix  B  is 
an  explanation  of  the  codes  used  in  Appendix  C. 

Plotting  routines  were  used  to  display  property  profiles  for  each 
station:   temperature,  salinity,  sound  speed,  and  density  (a  ) .   These  are 
compactly  plotted  four  stations  per  page  and  displayed  in  Appendix  D.   Sta- 
tions taken  in  the  deep  water  of  the  Barrow  Canyon  are  shown  two  per  page. 

IV.  RESULTS 

The  array  of  stations  occupied  is  shown  in  Figure  1  together  with  an 
ice-margin  position  based  principally  upon  observations  made  at  the  times 
stations  were  occupied.   The  ice-margin  is  thus  not  a  single  synoptic  view, 
but  a  progressively  distorted  one  which  is  more  useful  in  describing  ice- 
related  phenomena.   Synoptic  views  are  also  available. 
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An  interesting  feature  of  the  ice  in  the  region  of  Pt.  Barrow  is  the 
large  embayment  to  the  WNW.   Most  of  the  satellite-derived  observations 
show  the  embayment  more  or  less  closed  near  Pt.  Franklin  but  the  ship  found 
an  open  passage.   A  projection  of  ice  toward  Pt.  Franklin  is  a  feature  which 
persisted  for  several  weeks. 

Two  groups  of  measurements  were  made  from  a  hovering  helicopter,  one 
extending  about  kO   nmi  into  the  ice  and  another  group  of  stations  in  more 
or  less  open  water  which  supplemented  the  ship-based  stations.   Helicopter 
operations  proved  highly  useful,  but  could  not  always  be  used  where  they 
were  most  desirable  because  of  persistent  poor  visibility. 

The  distribution  of  f inestructure  is  shown  in  Figure  2  together  with 
isotherms  of  the  maximum  temperature  in  the  water  column.   The  intensity 
of  the  f inestructure,  measured  from  the  maximum  peak-to-peak  temperature 
fluctuation  in  the  profile,  is  categorized  as  weak,  moderate  or  strong. 
Weak  structure,  with  amplitude  0.2°  to  0.5°C,  is  indicated  by  a  circle  with 
a  central  dot;  medium  structure,  with  amplitude  0.5°  to  1.0°C,  is  indicated 
by  a  circle  with  a  cross;  strong  structure,  with  amplitude  greater  than  1°C, 
is  indicated  by  a  solid  circle.   Finestructure,  so  classified,  is  construed 
as  being  in  or  below  the  thermocline,  not  in  the  usually  shallow  upper  layer. 
An  open  circle  indicates  fluctuations  of  microstructure  size  or  the  absence 
of  any  notable  structure  at  all. 

We  decided  to  include  the  "nose"  classification  this  year,  after  the 
opposite  decision  in  1975,  because  of  a  feeling  that  it  may  be  an  inter- 
mediate step  in  forming  some  types  of  deeper  finestructure.   The  phenomenon 
we  have  called  a  "nose"  in  the  past,  in  its  simplest  form,  is  a  near-surface 
temperature  maximum  presumably  formed  by  the  cooling  of  an  upper,  warm 
mixed  layer  by  ice  or  the  overriding  of  such  a  layer  with  cold,  dilute  ice- 
melt  water  from  nearby  ice.   The  noses  this  year  differed  from  most  past 
results  in  that  they  were  sometimes  as  much  as  30  meters  thick  instead  of 
the  customary  10-15  m.   The  nose  may  be  complex,  having  various  levels  of 
finestructure  within  itself.   This  finestructure  sometimes  cuts  into  the 
usually  sharp  thermocline  so  as  to  make  the  depth  of  the  thermocline  doubt- 
ful and  the  nose  can  thus  grade  into  a  situation  which  we  should  classify 
as  normal  finestructure.   Further,  it  seemed  illogical  to  insist  that  the 
cooling  be  only  at  the  surface  since  an  initially  cooled  surface  skin  could 
be  subsequently  heated  by  solar  radiation  or  overriden  by  warm  dilute  water. 
Therefore,  a  nose  is  recognized 

a.  if  there  is  a  well-defined  surface  layer  5  or  more  meters  thick 
with  a  cooling  at  the  surface  of  at  least  0.5°C  as  compared  to 
the  maximum  temperature  in  the  layer. 

b.  if  there  is  a  well-defined  surface  layer  in  which  there  is 
moderate  or  strong  finestructure. 

The  nose  is  indicated  in  Figure  2  by  a  tab  on  either  side  of  the  circle. 
If  there  is  no  structure  in  the  nose,  the  tab  is  open;  if  there  is  medium 
or  strong  structure  in  the  nose,  the  tab  is  filled  solid. 
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It  will  be  noted  that  much  of  the  f ! nestructure  is  concentrated  be- 
tween the  2°  and  k°C    isotherms  of  maximum  temperature  in  the  water  column, 
a  conclusion  reached  in  studying  earlier  data.   However,  this  is  far  from 
being  an  absolute  rule:   there  are  regions  with  no  structure  between  the 
2°  and  *»°C  isotherms  and  one  instance  of  strong  structure  near  the  9° 
isotherm.   This  conclusion  therefore,  requires  closer  examination. 

The  following  general  results  are  presently  apparent. 

1.  The  ice  margin  crossings  in  the  vicinity  of  167°W  were  surprisingly 
like  the  corresponding  measurements  of  MIZPAC  72.   This  appears  to  be  a 
region  with  little  f inestructure. 

2.  Much  f inestructure  was  found  in  an  embayment  in  the  ice  roughly 
40  nmi  square  WNW  of  Barrow,  even  under  ice.   This  is  in  agreement  with 
the  findings  of  MIZPAC  71  which  operated  in  much  the  same  area,  but  with  a 
poorer  understanding  of  the  ice  distribution.   In  the  Barrow  Canyon,  where 
the  water  is  deeper,  f inestructure  was  found  as  deep  as  50  m.   This  is  15 
to  20  m  deeper  than  found  on  previous  cruises.   The  f inestructure  and  the 
embayment  probably  were  due  to  the  warm  water  of  the  coastal  current  which, 
in  this  case,  had  spread  westward  near  Pt.  Barrow.   It  is  interesting  that 
the  ice  was  fairly  close  to  shore  farther  south,  near  Pt.  Franklin,  and 

not  melted  away  as  one  might  expect  if  a  continuous  warm  current  moved  north- 
ward near  the  surface. 

3.  The  coastal  current  was  warmer  than  in  the  past,  up  to  8°  to  10°C. 
It  was  close  to  the  coast  between  Pt.  Franklin  and  Pt.  Barrow,  but  farther 
away  farther  south.   This  reinforces  our  feeling  that  the  tendency  of  the 
current  to  be  away  from  shore  in  the  more  southerly  latitudes  is  a  geo- 
graphical rather  than  a  temporal  phenomenon.   Five  crossings  of  the  coastal 
current,  between  Icy  Cape  and  Pt.  Barrow,  were  made  during  this  cruise. 
Analysis  of  these  crossings  should  provide  insight  into  the  route  of  the 
warm  Alaskan  Coastal  Current  as  it  flows  northeastward  from  the  shallow 
waters  off  Icy  Cape  to  the  deeper  waters  of  Barrow  Canyon. 

h.      The  salinity  of  bottom  water  behind  the  ice  edge  was  relatively 
high,  33A-33.6°/00,    like  in  MIZPAC  71  and  MIZPAC  72.   What  causes  this 
large  year-to-year  variation  in  mean  bottom  water  salinity,  often  greater 
than  0.5°/oo»has  not  been  established.   It  could  be  controlled  by  the  amount 
of  brine  rejection  occurring  in  winter  in  the  shallows  of  the  Chukchi  Sea, 
a  phenomenon  which  is  a  function  of  the  severity  of  the  winter.   However, 
the  cause  might  also  be  in  river  flows  into  the  Bering  Sea  during  the  pre- 
vious summer. 

At  the  extreme  bottom  of  many  casts  it  was  noted  there  was  an  apparent 
increase  in  salinity  in  the  last  half  meter  above  bottom.  We  are  undecided 
about  whether  this  increase  is  an  instrument  artifact  or  a  real  phenomenon. 
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APPENDIX  A 
DESPI KING  AND  DATA  EDITING 


Introduction 

Salinity  spikes  appear  to  be  considerably  more  serious  in  1977  than 
they  were  in  1975.   This  must  be  due  to  the  instrument  modifications  made 
in  the  interim.   Some  objective  despiking  routine  obviously  was  required. 
In  October  197**  we  had  tested  a  despiking  program  for  an  old  model  Plessey 
STD  based  entirely  on  the  concept  of  a  first-order  response  equation  for 
temperature  and  some  phase  shift  between  salinity  and  temperature,  which 
Scarlet  reported  using,  without  the  phase  shi ft,  in  1975.   The  method  did 
not  work  well  on  our  data.   Constants  which  would  correct  the  top  half  of 
a  spike  would  overcorrect  the  bottom  half.   The  reasons  may  lie  in  second- 
order  effects,  perhaps  complicated  by  the  fact  that  we  worked  from  the 
graphically  recorded  output  rather  than  from  digital  data  and  thus  included 
recorder  response  in  the  data.   Or  it  may  be  that  the  technique  which  Scarlet 
found  suitable  for  small  spikes  could  not  handle  the  very  large  spikes 
occurring  in  the  Chukchi  Sea  data  which  were  further  exaggerated  by  the  long 
thermometer  time  constant  of  the  STD.   Nevertheless,  necessity  prompted  us 
to  try  again,  with  some  success. 

The  data  editing  program  to  be  described  also  is  not  completely  new. 
We  had  one  in  MIZPAC  75.   But  the  present  program  is  considerably  more 
versatile  and  automatic.   It  contains  a  "ratchet"  or  "latch"  subroutine 
which  prevents  depth  reversals  similar  to  the  one  described  by  Scarlet. 
Our  routine  differs  in  that  Scarlet  threw  out  offending  data  whereas  we 
replace  offending  points  with  interpolated  points.   This  is  necessary  be- 
cause we  do  not  record  time  on  the  tape  and  must  assume  that  the  time  step 
is  uniform  when  despiking. 

Despiking  experiments  have  proceeded  at  the  Applied  Physics  Laboratory, 
University  of  Washington  (APL)  and  at  the  Naval  Postgraduate  School  (NPS) 
with  considerable  general  communication  between  the  two  laboratories.   The 
method  of  handling  the  short  time  constant  below  is  due  mainly  to  G. 
Garrison  of  APL,  the  method  of  correcting  temperature  and  the  long  time 
constant  is  due  to  R.  G.  Paquette  of  NPS.   The  details  of  implementation 
by  APL  and  NPS  are  different  and  we  have  not  intercompared  results.   The 
editing  program  described  was  evolved  at  NPS.   APL  uses  a  different  technique 

Facts  about  the  Instrument 

The  CTD  is  a  portable  digital  instrument  with  a  conventional  three- 
electrode  conductivity  cell  and  a  thermistor  as  a  temperature  sensor.   The 
cell  constant  of  the  cell  is  remarkably  stable,  judged  by  classical  ideas 
about  polarization  phenomena  in  cells  with  unplatinized  electrodes  and 
relatively  high  electrode  current  densities.   Perhaps  this  is  due  to  the 
high  operation  frequency,  ca_  10  kHz,  and  the  low  electrode  voltage,  13  mV. 
The  three  sensors  of  the  instrument  produce  electrical  frequencies  which 
are  counted  in  the  instrument  and  the  counts  are  recorded  digitally  on 


cassette  tape.   The  sensor  data  set  is  sampled  about  3  times  per  second 
in  the  order  conductivity,  temperature,  depth.   The  fact  that  the  conduc- 
tivity and  temperature  measurements  are  separated  by  about  30  ms  is  a  fac- 
tor in  the  generation  of  salinity  spikes  and  their  subsequent  removal. 

Salinity  Spikes-General 

Spurious  salinity  spikes  in  a  CTD  record  result  from  the  lack  of 
simultaneity  between  the  conductivity  and  temperature  measurement.   The 
spikes  are  of  consequence  only  when  there  are  sharp  temperature  transients. 
A  time  lag  may  occur  because  of  sampling  lag,  as  above,  or  it  may  occur  be- 
cause of  response  lags  in  the  sensors.   The  effect  exists  whether  the  salinity 
is  computed  digitally  or  by  an  analog  circuit.   When  the  temperature  is 
retarded  compared  to  the  conductivity  by  either  mechanism,  the  spike  has 
the  opposite  sense  from  the  first  derivative  of  the  temperature-time  curve 
and  conversely.   If  both  measurements  are  effectively  simultaneous  but  both 
are   wrong  because  of  sensor  lag,  the  salinity,  if  slowly  changing,  may  be 
correct. 

Sources  of  spikes 

In  our  investigation  of  the  cause  of  salinity  spikes,  which  sometimes 
exceeded  l°/0o  in  magnitude,  we  found  several  sensor  response  errors. 

°  A  first-order  lag  in  the  temperature,  with  time  constant  0.05  sec. 

°  A  first-order  lag  in  the  flushing  of  the  conductivity  cell  with  length 
constant  about  18  cm. 

°  An  error  due  to  non-simultaneous  sampling  of  conductivity  and  temperature 
and  to  a  physical  vertical  displacement  between  the  conductivity  cell 
and  the  temperature  sensor. 

°  A  long  time  constant  in  conductivity  response  apparently  due  to  heat 
storage  somewhere  in  the  cell  structure. 

Recognition  and  Evaluation  of  the  Long  Time  Constant 

It  may  be  of  interest  to  know  how  the  long  time  constant  came  to  be 
recognized  and  evaluated.   It  was  found  only  because  data  were  recorded  on 
the  upward  traverse  of  the  CTD  as  well  as  when  it  was  going  down.   It  was 
then  noted  that  the  down-going  salinity  curve  often  did  not  agree  with  the 
up-going  curve  in  the  region  below  the  thermocline  where  the  temperature 
was  unusually  constant.   No  temperature-induced  salinity  error  should  have 
been  generated  in  si  tu.   Obviously  the  error  must  have  been  a  long-lasting 
effect  of  a  salinity  or  temperature  transient  which  occurred  at  a  shallower 
level  and  affected  the  down  trace.   The  upward-going  trace  must  have  been 
essentially  correct  (except  for  the  neglected  effect  of  salinity  change  on 
conductivity)  so  long  as  the  temperature  was  constant  because  the  sensors 
usually  were  held  near  bottom  for  perhaps  10  seconds  before  beginning  to 
hoist.   The  phenomenon  under  discussion  may  be  seen  in  Figure  3,  which  shows 
the  down  and  up  salinity  traces  in  Station  80  as  well  as  the  down  and  up 
temperature  traces  and  the  corrected  salinity  trace.   One  may  easily  find  a 
time  constant  of  about  5  sec  in  the  approach  of  the  salinity  down  trace  to 
the  up  trace  below  the  thermocline,  taking  the  lowering  rate  to  be  1  m/sec. 
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If  such  a  time  constant  exists  in  other  STD's  and  CTD's,  there  must  be 
difficulties  in  removing  spikes  by  the  use  of  a  single  time  constant.   This 
is  because  the  short  and  long  time  constants  contribute  differently  to  narrow 
spikes  than  to  broad  spikes;  the  short  time  constant  influences  the  former 
more  than  it  does  the  latter  and  conversely  for  the  long  time  constant. 

The  Possibility  of  Second-Order  Effects 

All  of  the  discussion  of  spike  correction  below  will  assume  the  simple 
first-order  response  equation  for  sensors.   It  should  be  kept  in  mind, 
however,  that  substantial  second-order  components  may  be  present.   Whenever 
there  are  two  thermal  masses  in  series  with  two  thermal  resistances  involved 
in  the  transfer  of  heat  from  the  water  to  the  interior  of  the  sensor,  a 
second-order  term  enters  the  response  equation.   The  temperature  and  salinity 
sensors  must  have  had  these  characteristics  to  some  extent,  but  it  was  not 
possible  to  find  the  multiplier  of  the  second-order  term.   A  second-order 
corrector  basically  takes  out  curvatures.   It  was  found  that  the  second-order 
constant  could  be  adjusted  to  take  out  anomalies  in  the  curves  beautifully, 
but  there  was  no  way  to  tell  whether  the  corrections  corresponded  to  reality. 
Therefore,  after  a  few  experiments,  the  refinement  was  abandoned.   It  nay 
be  possible  to  revive  it  at  some  later  date  when  more  is  known  about  the 
causes  of  variability  of  the  sensor  response  constants. 

Lack  of  Constancy  of  the  "Constants" 

A  major  problem  in  devising  a  spike  correction  technique  is  the  fact 
that  the  response  constants  of  the  sensors  are  not  constant.   They  change 
with  rate  of  lowering,  wire  angle  and  relative  horizontal  motion  between  the 
sensor  head  and  the  surrounding  water.   If  the  constants  had  been  nicely 
fixed,  there  would  have  been  straightforward  techniques  for  evaluating 
several  constants.   As  it  was,  there  was  little  hope  of  evaluating  the 
second-order  constant  and  little  justification  for  trying  to  compute  the 
effect  of  the  flushing  constant  of  the  cell  upon  the  conductivity  directly, 
a  natter  which  will  be  discussed  below  when  the  methods  of  spike  correction 
are  related  in  detail. 

Theory  and  Algorithms  for  Correction  of  Spikes  -  Thermometer  Correction 

In  the  following  paragraphs  the  nature,  causes  and  method  of  correction 
of  the  several  tine  constants  are  discussed. 

We  assumed  that  the  thermometer  followed  a  sinple  first-order  law 

t  -  r  -  k   £1 

'         Tdt 

where  T  is  the  water  temperature,  T'  the  observed  temperature,  k   the  tine 
constant  of  the  thermistor  and  t_  the  time.   The  bare  thermistor  had  a  noninal 
time  constant  of  0.03  sec.   The  protective  shroud  was  expected  to  add  to  this 
After  a  number  of  trials  we  picked  0.05  sec  as  the  largest  value  which 
would  not  produce  an  overshoot  at  the  bottom  of  the  sharpest  therr->ocl  i  nes . 
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The  temperature  derivative  was  obtained  from  the  first  central  difference 


/  dTM    T'  .  -  Ti 


where  h^  is  the  time  step  between  samples,  0.33  sec. 

The  Time  Lags  and  the  Flushing  Time  Constant 

The  difference  in  sampling  times  of  temperature  and  conductivity  was 
0.03  sec  (temperature  latest)  plus  an  addendum  of  0.01  sec  because  the 
thermometer  was  1  cm  deeper  than  the  mouth  of  the  cell.   These  were  com- 
pensated by  interpolating  backward  toward  the  previous  temperature  a  frac- 
tional distance  (0.01  +  0.03)/h  or  0.012,  using  the  equation 

AT.  =  k  T.  .  +  (1-k  )  T.  -  Ti 
j    c  j-1      c   j    j 

where  the  constant  k  would  be  expected  to  be  0.12,  from  above,  but  was 
more  effectively  about  0.65.   The  extra  0.53  can  be  shown  to  approximate 
a  correction  for  the  flushing  lag  in  the  conductivity  cell,  which  apparently 
has  a  time  constant  of  0.53  h  or  0.18  sec.   This  corresponds  to  a  length 
constant  of  18  cm,  or  one  cell  length,  which  is  reasonable. 

There  are  some  known  theoretical  weaknesses  in  treating  the  flushing 
time  constant  of  the  cell  in  this  way.   The  method  involves  the  assumption 
that  there  are  no  significantly  large  conductivity  gradients  that  are 
not  caused  exclusively  by  temperature  changes,  which  is  the  same  as  assuming 
that  the  true  salinity  always  changes  slowly.   It  also  involves  the  assumption 
that  the  first  backward  temperature  difference  is  proportional  to  9c/8t, 
where  c_   is  the  conductivity.   These  are  not  always  good  assumptions.   However, 
reasonable  success  with  the  techniques  led  us  to  postpone  further  refinements 
to  a  later  date. 

The  Long  Time  Constant 

The  long  time  constant  in  the  cell  was  modeled  successfully  as  a  thermal 
mass  coupled  to  true  water  temperature  by  a  thermal  resistance  corresponding 
to  a  time  constant  k.  .   This  mass  passes  fraction  f_  of  the  difference  be- 
tween its  temperature  and  true  water  temperature  to  the  cell.   It  is  a 
fairly  realistic  model.   The  thermal  mass  might  be  the  cell  body  itself  or 
the  surrounding  protective  shroud. 


The  temperature,  TC,  of  the  thermal  mass  is  accumulated  at  each  step, 
starting  with  TC.  =  T.  and  assuming  that  the  mean  temperature  difference  be- 
tween TC  and  T  drives  the  temperature  change  in  accordance  with 


TC.  =  TC.  .  + 
J     J-1 


T.+T.  . 
J   J'1 


TC.+TC.  . 

_J r± 


which  ultimately  gives 
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TC.  = 
J 


T.  .  +  T..+  TC. .  .  (2kL  -  h" 
J-l     J      J"'   u 


2KL+h 


Then  the  temperature  addendum  contributed  to  the  cell  is  ATC .  =  (TC.-T.)F, 
where  F  is  about  0.15-  J      J   J 

Difference  Equation  for  Converting  Temperature  Error  to  Salinity  Correction 

At  this  point  in  the  discussion,  all  the  corrections  have  been  com- 
puted in  terms  of  their  equivalent  temperature  corrections  assuming  C.  may 
be  taken  as  the  observed  conductivity.   The  effects  of  k  and  k  on  salinity 
were  then  computed  separately  by  a  difference  equation  derived  from  the 
salinity  temperature  -  conductivity  tables  of  H.O.  Pub.  SP-68  in  the  vicinity 
of  2°C  and  30°/oo.   This  equation  is 


S-S'  = 


-11.0  X  10~5  (T-T1)  -  2.28  X  10~5S'(T-T') 
76.8  X  10~5  +  2.28  X  10"5  T 


When  the  short  time  constant  and  the  time  lag  were  being  corrected,  T-T1 

was  replaced  by  AT.  above.   When  the  long  time  constant  was  being  corrected, 

T-T'  was  replaced  Ay  T.  +  ATC.  -  T'.. 

J      J    J 

Example  of  Corrections 

Figure  k   shows  the  original  salinity  curve  from  Station  80,  the  two 
components  of  the  despiking  correction  and  the  resultant  corrected  curve. 
There  were  nearly  always  some  small  irregularities  left  after  despiking. 
These  were  smoothed  by  means  of  a  9~point  running  mean.   The  corrected 
salinity  curve  before  smoothing  may  be  seen  as  the  dotted  curve  in  Figure  3- 

Method  of  Finding  Despiking  Constants 

The  constants  for  despiking  were  determined  in  the  following  way.   The 
despiking  program  was  used  intensively  on  about  ten  different  stations, 
including  some  of  the  up  traces,  to  become  familiar  with  the  effects  of  the 
several  constants.   In  this  way  the  general  range  of  the  constants  and 
their  variability  was  found.   The  constants  k.  and  F  were  determined  as 
those  which  would  make  the  down  trace  and  the  up  trace  match,  or  be  nearly 
parallel  below  the  thermocline.   The  constant  k  was  chosen  to  remove  spikes 
as  well  as  possible.   Sometimes  it  was  necessary  to  adjust  F  also  to  remove 
spikes  which  meant  that  perfect  matching  of  up  and  down  traces  was  not 
possible.   However,  perfect  matching  is  not  a  requirement  unless  the  salinity 
curve  is  nearly  vertical  near  bottom  because  a  salinity  gradient  causes  a 
small  error  in  cell  response  which  is  not  corrected  by  k  ,  as  was  explained 
previously.   The  despiking  program  was  then  incorporated  into  the  general 
editing  program.   During  the  general  editing,  further  adjustments  of  the 
despiking  constants  were  made  before  writing  the  corrected  values  on  a  new 
tape.   A  description  of  the  general  editing  program  follows. 

Some  Tests  of  Responses  at  Sea 

Two  stationswere  duplicated.   In  the  first  run  of  Station  46,  the  lower  ing 
was  made  by  backing  up  1/2  meter  for  each  meter  of  lowering  through  half  of 
the  water  column.   The  rest  of  the  drop  was  made  at  normal  speed,  ca  lm/sec. 
In  the  second  run  the  lowering  was  normal.   It  was  expe.cted  that  the  response 
of  the  sensors  would  be  notably  different.   The  spiking  was  less  serious 
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Figure  4.  Station  80.  The  down-going  salinity  profile  is  on  the  right  as 
a  solid  curve.  The  dashed  curve  is  the  result  of  despiking  and  smoothing. 
The  two  curves  to  the  left  give  the  contributions  of  the  short  time  constant 
and  long  time  constant  corrections  to  the  despiking.  See  Figure  3  for  the 
unsmoothed  corrected  curve  and  to  see  the  temperature  transients  which 
caused  the  salinity  spikes. 
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In  the  first  case.   However,  it  is  curious  that  the  despiking  constants 
required  for  the  first  were  the  normal  ones,  k  =0.65,  F=0.15  whereas,  for 
the  second,  k  =0.38  and  k  =0.115  were  required.   The  corrected  results  are 
very  similar  except  near  the  surface  where  rapid  real  changes  in  properties 
are  common.   In  Station  107,  the  first  drop  is  normal  and  the  second  was 
made  with  the  thermistor  taped  onto  the  multiplexer  housing  facing  upward. 
The  smaller  spikes  of  temperature  and  salinity  were  smoothed  in  the  second 
case,  presumably  by  the  turbulence  generated  by  the  multiplexer  housing. 
This  is  a  station  in  which  the  temperature  transients  are  moderate  and, 
although  there  were  differences  in  the  original  curves,  they  are  not  dramatic 
enough  to  interpret  in  terms  of  the  different  cell-to-thermistor  distance. 
Both  lowerings  were  corrected  with  the  standard  constants  and,  when  corrected, 
looked  very  similar. 

Requirements  for  General  Editing  of  Data 

The  CTD  data  require  considerable  editing  besides  the  despiking,  as 
will  be  noted  below. 

°  Spurious  high  salinities  frequently  occur  near  the  surface  on  the 
down  trace,  and  the  first  one  or  more  depths  are  occasionally  very 
large.   High  apparent  salinities  are  created  by  very  low  actual  con- 
ductivities due  either  to  melt  water  or  to  bubbles  in  the  cell.   The 
system  records  only  four  digits  of  a  five-digit  frequency  count. 
Therefore,  foldover  occurs  as  a  result  of  low  values  and  the  last 
four  digits  are  at  the  top  of  the  range  and  are    read  as  high  values. 

°   Occasionally  an  absurd  single  value  may  occur  in  any  of  the  measured 
parameters.   Sometimes  the  error  is  too  small  to  be  distinguished 
automatically  from  real  fluctuations. 

°   Near  the  top  of  a  lowering  there  are  duplicate  values  generated 
during  the  "soaking"  period  and  also  while  the  sensors  are   waiting 
near  bottom  while  buttons  are  pressed  on  deck.   The  shallowest  depth 
may  be  interspersed  among  other  depths  near  surface;  a  similar  result 
occurs  near  bottom  because,  after  striking  bottom,  the  sensor  was 
immediately  raised  one  or  two  meters.   A  method  is  required  to  re- 
move the  unwanted  records  and  to  start  as  shallow  as  possible  and 
end  as  deep  as  possible. 

°   Some  lowerings  were  so  full  of  noise  as  to  require  replacing  them 
with  the  corresponding  up  trace.   In  such  cases,  special  adaptations 
of  the  noise-removal  and  despiking  routines  were  necessary  and  the 
data  set  was  then  inverted  top-to-bottom. 

°  When  the  ship  rolls,  loops  are  generated  in  the  recorded  traces 
because  of  sensor  response  problems.  Records  distorted  by  such 
reversals  in  depth  had  to  be  adjusted  before  despiking. 

°  The  several  sensors  have  calibration  errors,  additive  in  the  case 
of  temperature  and  salinity,  and  a  factor  for  converting  from  the 
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fresh-water  to  a  salt-water  calibration  in  the  case  of  depth.   A 
substantial  depth  discrepancy  will  be  noted  between  down  and  up 
curves  in  Figure  3.   The  cause  is  only  partly  known  and  it  was  not 
corrected.   The  major  contributions  are  a  sampling  lag  of  twice 
255  ms  and  the  probability  that  temperature  features  in  the  water 
were  carried  upward  by  the  multiplexer  housing  of  the  CTD  which 
preceded  the  sensors  on  the  way  up. 

Description  of  the  General  Editing  Program 

The  complete  data  editing  program  is  a  group  of  routines  which  may  be 
invoked  in  sequence  by  commands  entered  on  a  control  card  read  immediately 
before  the  data  to  be  edited  were  read  from  tape.  The  following  functions 
could  be  performed. 

1)  Eliminate  an  entire  "station".   Here  the  word  station  means  those 
data  isolated  on  the  tape  by  interrecord  gaps.   Usually  it  was  the 

up  trace  which  was  to  be  eliminated. 

2)  Eliminate  sequential  data  records  in  up  to  two  places  on  any 
station,  the  places  designated  by  beginning  and  ending  serial  numbers 
of  the  records.   This  served  to  remove  faulty  values  frequently  found 
at  the  beginning  of  a  lowering  and  the  repetitious  values  usually 
found  where  the  sensor  head  was  stopped  on  or  near  bottom. 

3)  Interpolate  between  two  good  records  as  many  records  as  were 
previously  present  in  a  faulty  intervening  group.  This  could  be 
done  in  two  places  on  any  station. 

k)      Replace  up  to  three  records  with  images  punched  on  cards. 

5)  Remove  single-point  spikes  in  depth,  temperature  or  salinity.   A 
single-point  spike  found  in  one  of  these  was  almost  certainly  an  arti- 
fact and  was  replaced  by  the  median  between  the  J-l-th  and  J+l-th 
measurement. 

6)  Apply  a  depth  ratchet  so  that,  after  the  first  20  points,  the 
depth  cannot  decrease.   Records  in  which  the  depth  has  decreased  as 
compared  to  D.  are  replaced  by  an  equal  number  of  records  interpolated 
evenly  between  D.  and  the  next  depth  which  is  equal  or  greater.   At 
the  beginning,  tne  routine  automatically  discards  records  prior  to 
the  minimum  depth  found  in  the  first  20  records.   Where  this  routine 
was  to  be  applied  to  an  up  trace,  a  variable,  UP,  was  set  true  which 
caused  all  the  depths  to  be  temporarily  replaced  by  their  negatives 

so  that  the  ratchet  routine  would  work. 

7)  Apply  the  despiking  routine.   The  constants  could  be  changed  for 
each  station. 

8)  Invert  the  sequence  of  records,  putting  the  top  at  the  bottom, 
this  for  a  few  cases  in  which  the  down  trace  was  too  faulty  to  be 
used. 

9)  Make  additive  corrections  to  salinity  and  temperature  based  upon 
the  comparisons  with  Nansen  bottles  and  correct  depth  for  water  density. 
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10)   Recompute  sound  velocity  and  sigma-t  from  the  corrected  salinity 
and  temperature. 

The  program  then  wrote  the  corrected  data  on  a  new  tape,  produced  a 
printer  plot  of  the  corrected  and  uncorrected  salinities  and  the  contribu- 
tions of  the  long-time  constant  and  lag  corrections  separately. 


17 


APPENDIX  B 
EXPLANATION  OF  HEADING  CODES 


The  heading  of  the  printed  output  uses  the  coding  and  format  from 
NODC  Publication  M-2,  August  1964,  with  a  few  exceptions.   Heading  entries 
which  are  not  self-explanatory  are  as  follows:   MSQ  is  the  Marsden  Square, 
and  DPTH  is  the  water  depth  in  meters.   Wave  source  direction  is  in  tens 
of  degrees,  but  the  direction  99  indicates  no  observation.   The  significant 
wave  height  is  coded  by  Table  10  (Code  f  2  ~  height  in  meters)  and  the  wave 
period  is  coded  by  Table  11  (Code  t  2  ~  period  in  sec);   in  each  case  X 
indicates  no  observation.   Wind  speed,  V,  is  coded  as  Beaufort  force,  Table 
17.   The  barometer  is  in  millibars,  less  1000  if  more  than  3  digits;  wet 
and  dry  bulb  temperature  are  in  degrees  C.   The  present  weather  is  from 
Table  21  with  cloud  type  and  amount  from  Tables  25  and  26,  respectively. 
The  combination  4X9  indicates  that  clouds  cannot  be  observed  usually  be- 
cause of  fog  conditions.   The  visibility  is  from  Table  27  which  is  roughly 
in  powers  of  two  with  Code  4=1-2  km.   The  ice  concentration,  IC,  is  in 
oktas;  amounts  less  than  1  okta  are  preceded  by  a  minus  sign  and  indicate 
concentrations  in  powers  of  ten,  e.g.,  10   =  -4. 

The  entry,  COD,  is  a  code  to  identify  the  accuracy  of  each  station 
position  based  upon  the  navigation  system  used.   Code  1  indicates  a  posi- 
tion determined  by  visual  sightings  or  radar,  Code  2  a  position  determined 
by  navigation  satellite,  and  Code  3  a  position  determined  by  DR. 
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APPENDIX  C 
HEADING  DATA  FOR  MIZPAC  77  STATIONS 


Heading  data  are  listed  on  the  following  pages  for  MIZPAC  77.   The 
coding  conventions  are  those  described  in  Appendix  B.   Note  that  Stations 
67H  through  7 1 H  and  97  are  missing.   Other  stations  in  the  helicopter 
series  have  much  of  the  heading  information  missing. 
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APPENDIX  D 
PROPERTY  PROFILES  FOR  MIZPAC  77  STATIONS 

This  section  contains  plots  of  temperature,  salinity,  sound  velocity 
and  sigma-t  for  all  of  the  stations  of  MIZPAC  77  which  were  recovered 
from  the  cassette  tapes  successfully.   Station  97  is  missing  and  of  the 
sequence  65H  through  71 H ,  allotted  for  a  helicopter  expedition,  only  the 
first  two  were  successful.   Other  "H"  suffixes  also  indicate  helicopter 
stations.   It  will  be  noted  that  Stations  h6   and  107  are  duplicated.   Ex- 
planations of  the  differences  between  duplicates  are  given  in  Appendix  A. 

The  basic  four-per-page  plot  has  a  maximum  depth  of  70  m.   All  the 
stations  were  plotted  in  this  way.   In  addition,  deeper  stations  were  plot- 
ted on  a  140-meter  depth  scale,  two  per  page.   These  are  interleaved  with 
the  smaller  plots.   To  assist  in  distinguishing  curves  the  salinity  curve 
has  been  doubled  with  a  line  0.01  inches  to  the  left  and  the  temperature 
has  been  treated  similarly  with  0.015  inch  spacing.   The  curves  are  also 
labeled,  T  for  temperature,  S  for  salinity,  SV  for  sound  velocity  and  ST 
for  sigma-t. 
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